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(Dated: July 24, 2006) 

Abstract 

We present preliminary measurements of D ^ -^l^ ^-^id D Ks'^ branching fractions using 
281 pb-i of *" (3770) data at the CLEO-c experiment. We find that H(D° ^ K^tt^) is larger 
than B{D^ Kl'^^), with an asymmetry of R{D°) = 0.122 ± 0.024 ± 0.030. For B{D+ K^tt^) 
and B{D^ i^^vr"'"), we observe no measureablc difference; the asymmetry is R{D~^) = 0.030 it 
0.023 lb 0.025. Under reasonable theoretical assumptions, these measurements imply a value for 
the — > K^TT^ strong phase that is consistent with zero. The results presented in this document 
are preliminary. 



*Submitted to the 33'''^ International Conference on High Energy Physics, July 26 - August 2, 2006, Moscow 
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I. INTRODUCTION 



The meson, with quark composition cd, decays by a Cabibbo-allowed decay to K^n^ 
{dc dsWy {sd){ud) =^ K^n^), and by a doubly-suppressed decay to i^°7r+ {dc — *• 
ddWy —>■ ddus {ds){ud) =^ K^tt^). Similarly, the meson, with quark composition cu, 
decays by a Cabibbo-allowed decay to K'^tt^, and by a doubly-suppressed decay to K'^n^. 
The observable final states are not or K^, but rather or K^. As pointed out by 
Bigi and Yamamoto many years ago 0, interference between Cabibbo-allowed and doubly- 
suppressed transitions leads to differences in the rates for D —* K^tt and D Kgir. Here we 
describe a search for differences in the decay rates, both for Kgir^ vs. K^-k^ 

and for K^ti^ vs. K^n^. (Throughout, charge-conjugate modes implied, except 

where noted.) 

For these measurements we have used a sample of 281 pb^^ e^e^ ^"(3770) events, 
produced with the CESR-c storage ring and observed with the CLEO-c detector. 

The data sample contains approximately 820,000 D~^D~^ events and 1,030,000 D^D^ 
events, as well as e~^e~ — > uu, dd, ss continuum events, e^e^ — > r"'"r~ events, Bhabha 
events, and /x-pairs. The resonance (3770) is below the threshold for DDtt, and so the 
events of interest, e^e^ ^" DD, have D mesons with energy equal to the beam energy 
and a unique momentum. 

For the decays D — »• K^tt, we make no attempt to detect the as this is not feasible 
with the CLEO-c detector. Rather, we fully reconstruct a tag D on "the other side," detect 
the 71, and compute the missing mass squared. Our signal is a peak at the mass squared. 

Explicitly, for D^ —>■ K^n^, we reconstruct the D~ in 6 decay modes. We do this by 
requiring that the candidate D~ has energy consistent with the beam energy, and "beam- 



constrained- mass" (y £'|g^^ — |SPjp) consistent with the D~ mass. Given a reconstructed 
D~ meson, we require that the remainder of the event (the "D"^ side") contain only one 
charged track (for the tt"*"), and that any calorimeter clusters do not form a 7r°. We then 
compute the missing mass Mx in the reaction e~^e^ D^tt^X. The result is shown in 
Figure n where the peak is evident. 

For the decay D^ K^t^^, we follow a similar procedure. We reconstruct D^ in the 
decay modes D'^ K^tt^ , K~^7r~7r^, and i^"'"7r~7r"'"7r~. Given a reconstructed D'^ meson, we 
require that the "D° side" contains no charged tracks, only one 7r° —>■ 77, and no rj 77. 
We then compute the missing mass Mx in the reaction e~^e~ —>■ D^tt^X. The result is shown 
in Figure 121 with the K^^ peak evident. 

For D~^ K^TT~^, we use the result of an independent CLEO-c analysis that measures 
many D° and D^ hadronic branching fractions, including D^ K'^n'^. This decay is 
directly reconstructed, using Kg tt+tt^. Unfortunately, the analysis on the full 281 pb~^ 
sample is not yet complete; the result we use is based on a 56 pb~^ subset. An updated 
result, with approximately half the uncertainty, will be available soon. 

No previous CLEO-c analysis has measured D° Kgir^, so we have also analyzed this 
mode. We directly reconstruct this mode, either tagged with a reconstructed D^ ("double 
tag") or not tagged ("single tag"). 
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II. QUANTUM CORRELATION IN DECAYS 

The situation for — > K^ti^ vs. K^ti^ has an added feature. When and D'^ are 
pair-produced through a virtual photon(J'^'^ = 1 ), they are in a quantum coherent state. 
Then the decays of and will follow a set of selection rules. They cannot decay to 
CP eigenstates with the same CP eigenvalue, if we ignore CP violation in D^/D^ system. 
On the other hand, decays to CP eigenstates with opposite CP eigenvalues are enhanced. 
Similarly, all other final states are subject to such inteference effects. As a result, the 
measured branching fractions in this / system differ from those of isolated mesons. 
The measured branching fractions of the same mode by double tag and single tag methods 
will also differ from each other, especially for CP eigenstate modes. 

The quantum correlation effects are shown in Table lU where " f " stands for flavored 
modes, " X " stands for everything, " 5"+ " stands for CP even modes, " S_ " stands for CP 
odd modes, and 

_ mi - m2 _ Pi - Pa {f\D^) _ 

^" P ^" 2P ''^'^ 

where is the amplitude ratio of "wrong sign" decay(D° K^tt~ for example) to "right 
sign" decay (D*^ K~ti^ for example) and 5 is the phase difference. 



TABLE I: Single Tag and Double Tag yields for C = — 1 D^D^ events, to leading order in x and y. 





s+ 


S- 


f 


NBfBs+{l + r) + 2rfCos6) 


NBfBs-il + rj- 2rfCos6) 


X 


2NBsJl-y) 


2NBs-{l + y) 



Since the CP eigenvalue of K^tt^ is odd and the CP eigenvalue of K^n^ is even we can 
see from Table H] that we will overestimate the K^tt^ branching fraction in the double tag 
method and underestimate the i^^vr" branching fraction. For single tag measurements, the 
effects are small since y is tiny. 

Our procedure is the following: 

1. We measure the "branching fraction" for -^5^°) untagged. This gives us 
B{D^ Kl-n^){l + y). Since y is very small, 0.008 ± 0.005[PDG], we can correct 
for it, obtaining K'^s^^). 

2. We measure the "branching fraction" for Kgir'^, with three different flavor tags. 
Each gives us ^s7r°)(l - 2rfCosS + rj). Using B{D^ ^s^°) obtained from 
the untagged measuement, we obtain 1 — 2rfCos6 + r'j, for each flavor tag. Since rj is 
known for these three flavor tags, from this we can compute 1 + 2rfCos6 + r'j for each 
flavor tag. 

3. We measure the "branching fraction" for — >■ K^ir^, with the same three flavor tags 
as those used for Each gives us B{D^ _ft'°7r°)(l + 2rfCosS + rj). Using 
the values of 1 + 2rfCos6 + rj for each tag obtained in (2), we obtain B{D^ — * K^tt^), 
for each of the three flavor tags. These three measurements are then averaged for the 
final result. 
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TABLE II: Yields in both Monte Carlo (top) and data (bottom). The Monte Carlo input branching 
fraction for K^tt^ is 1.06%. 



A/To Hp 


nO . z<'0„o 

s 


^signalregion 
i-XJ-j oiU-CUdliU. 


128407 
fifi9Q 7 




tJO ± .o 


*MHfir\5i'nri tin mI'tpi r'^'f^H 

kJlLlC UaillKJL o U. U Ui Ci\^ UCvJ. 




Nd + Nf) 


36295180 


Efficiency 


28.94% 


B 


1.053±0.007% 


^signalregion 

AE sideband 


8726 
944.4 


K'^ sideband 


294.3 


Sideband subtracted 


7487.2 




1013314 


B 


1.212±0.016% 



III. D K^TT MEASUREMENTS 

The value of B{D+ Ks^^) is taken from reference (1.55 ± 0.05 ± 0.06)%. 
B{D^ — > K^n^) is measured with two methods: single tag and double tag. 

A. Single Tag L>° ^ K'^tt^ 

Candidates for — > K^n^ in untagged events were formed by combining a Kg, recon- 
structed by a pair of charged tracks through the decay Kg —>■ n~^7i~ , and a 7r° from pairs of 
photons detected in the Csl crystal calorimeter. The invariant mass of Kg candidates was 
required to be within 3 standard deviations of the known Kg mass. The sideband of Kg mass 
is from 4 standard deviations to 7 standard deviations on both sides. The invariant mass of 
7r° candidates was required to be within 4 standard deviations of the known vr" mass. Pho- 
tons with energy less than 30 MeV were not considered. Both beam constrained mass and 
AE were required to be within 3 standard deviations of the nominal value. Two sideband 
subtractions were used to subtract the background. AE sideband subtraction was used to 
subtract the continuum and combinatoric background. Kg mass sideband subtraction was 
used to subtract the peaking background under AE and Mbc distributions. This peaking 
background was formed from a real D with the final state tt+tt^vt'^, in which 771(7? "•"vr") hap- 
pens to be within the Kg mass window. This mode is not decayed from a Kg resonant state, 
so the mass of tt+tt" shows a flat distribution. The Kg mass sideband subtraction removes 
this background effectively. 

All the yields of Monte Carlo and data in the signal region and the sideband region are 
shown in Table ITTl By using the luminosity and cross section of e'^e" D^D^ , we get the 
number of D's in data. Combining with the efficiency from signal Monte Carlo, we get the 
branching fraction for —>■ Kgir^ in the untagged data sample. 

The systematic uncertainties are listed in Table ITTTl The uncertainties due to vr*^ recon- 
struction efficiency will cancel in the comparison of branching fractions for -^lTt", 
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TABLE III: Systematic uncertainty for single tag 5(1)° 



AE cut 


0.5% 


iracKing ernciency 


U. 1 /o 


AE sideband 


0.82% 


Kg efficiency 


1.1% 


Kg sideband 


0.28% 


Cross section 


2.75% 


7r° efficiency 




All (without 7r° efficiency) 


3.20 



D° Kgir'^, we will not include that uncertainty here. The dominant uncertainty comes 
from the cross section, which is used to calculate the number of D's. This cross section is 
base on the 56 pb~^ dataset. Soon, when the 281 pb~^ result comes out, this uncertainty 
will improve. 

Combining all the results above, the single tag branching fraction for K^tc^, without 

7r° systematic uncertainty, is 1.212 ± 0.016 ± 0.039%. 

B. Double Tag i^^vr" 

For the tagged branching fraction of Kgir'^, was fully reconstructed as — >■ 

K+7r~, — > K^Tc^TT^, or K+tt^tt+tt^. In the tagged sample, we reconstructed 

— > KgTT^ in the same way as in the untagged case. All the requirements were unchanged, 
except there were additional requirements regarding the tag side. The tag Z)° was required 
to be within 3 standard deviations in both the AE and Mbc distributions. For tag mode 
— >• K~^TT~Tr^, the energy of the tag-side 7r*^'s lower-energy shower was required to be 
greater than 70 MeV. This requirement made the background in AE distribution flatter 
and thus more suitable for the use of AE sideband subtraction to get the number of D's in 
the tag side. 

In the tagged data sample, since the K^tt^ mode has relatively fewer neutral and charged 
particles than an average decay, it is easier to reconstruct the tag when — > K^tt^, 
especially for the tags with more charged or neutral particles. Therefore, the branching 
fraction of the signal mode is biased in the subset of the selected tag. By checking the 
Monte Carlo truth in the tagged sample, we obtained a correction factor for this tag bias. 

Tag side AE sideband subtraction was used to subtract fake D events. The signal side 
Kg mass sideband was used to subtract peaking background under signal side AE and Mbc 
distributions. Just as in the untagged case, the peaking background was formed from a real 
tag D with, on the signal side, the final state 7r~^7T~TT^, in which m(7r"*"7r~) happens to be 
within the Kg mass window. 

With the yields and effciencies from signal Monte Carlo, we computed the branching 
fraction in Table IIVI 

The uncertainties due to data and Monte Carlo differences in 7r° and Kg reconstruction 
efficiencies are the same as in the untagged case. These two uncertainties will cancel in the 
ratio of these two branching fractions. The uncertainties due to AE sideband subtraction 
and Kg sideband subtraction were also estimated in a similar way as in untagged case. 

The average of the results for the three tag modes is (1.032 ± 0.047)%, significantly 
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TABLE IV: Branching Fraction for /Cgvr" in double tag method, "s-s" means sideband- 

subtracted for tags, and background subtracted for signal yield, "a" is the correction factor for tag 
bias; Error is statistical only. The Monte Carlo input branching fraction for Kgir^ is 1.06%. 





Tags 


Signal Yield 


iVlu 


Corrections 




MC 


raw s-s 


raw s-s 


eff(%) 


a 


BR(%) 


Ktt 


859919 855772 


3119 2958.5 


33.04 


1.00 


1.046±0.020 


Kmr^ 


1186407 1147754 


4815 4129.5 


32.57 


1.014 


1.089±0.020 


KStt 


1268558 1228967 


4480 4156.5 


31.20 


1.033 


1.049±0.017 


Data 


raw s-s 


raw s-s 


eff(%) 


a 


BR(%) 


K-K 


48095 47440 


172 155 


33.04 


1.00 


0.989±0.088 


Ktt-k^ 


67576 63913 


248 203 


32.57 


1.014 


0.975±0.082 


KStt 


75113 71039.5 


276 256 


31.20 


1.033 


1.118±0.075 



different from the untagged result, illustrating the effect of quantum correlations. (Because 
the quantum correlation correction factor, 1 — 2rcosS + r^, is tag- mode-dependent, this 
average is not otherwise of interest.) 

IV. D kItt measurements 

We measure the D K^tt branching fractions with a missing mass technique. We 
reconstruct the tag Z) in 3 modes and 6 D~ modes, and we combine it with a 7r° or vr"*" 
to form a missing mass squared. The D K^tc signal is a peak at the K'^ mass squared 
((0.49772 GeV)2 = 0.24773 GeV^). 

To remove D K^ir events, as well as other backgrounds, we require that the event 
contain no extra tracks or 7r°'s beyond those used in the tag D and the vr. This veto removes 
about 90% of D — > KgH events and a few percent of D ^ Kin events. 

A. D+ KI-K+ 

We reconstruct tag D~'s in the 6 decay modes D~ —>■ K^tt~ti~ ^ K~^7i~7i~7t'^, K^tt^ , 
Kg7r~7r^, Kg7r~n~7r~^ , and K~^K~n~ . Candidates must pass AE and Mbc cuts. 

The tag reconstruction efficiency is generally higher when the signal decays to fT^vr"'" 
than for generic decays because Ki^tt~^ has only one charged particle and at most one 
calorimeter cluster. This biases the sample of tagged events in favor of signal events. There- 
fore, we include a factor in the branching fraction calculation to correct for this tag bias. 
The factor, measured in Monte Carlo, is the ratio of the tag reconstruction efficiency when 
the decays to i^^^vr"'" to the efficiency when it decays to anything else. 

The efficiency for observing K^n^, given that the tag was successfully recon- 

structed, is measured in signal Monte Carlo. It is essentially the efficiency for finding the 

71+ . 

The missing mass squared distribution, with all tag modes added together, is shown in 
Figure^ The lines show a fit to determine the signal yield; each line represents a background 
component added cumulatively. The most prominent feature is the signal peak at the 
mass squared (~0.25 GeV^). A number of backgrounds are also present. First, fake D~ 
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Missing mass squared (GeVA2) 



FIG. 1: Fit for ^l'^^ yield using all tag modes. The many colored lines represent the various 

background components, added cumulatively. The green peak is the contribution of ^5^^ 
events that are not removed by the extra track and vr^ vetoes. 

candidates produce a background which is estimated from an Mbc sideband. All of the other 
backgrounds come from other decays. The largest of these are Kgir'^ (green peak 

under the signal), rjiT'^ (peak on the right-side tail of the signal), 7r°7r+ and /i"*"//^ (peak on the 
left of the plot), K^7i~^7i^, and vr+vr^vr^. The shapes and efficiencies of these backgrounds are 
determined from Monte Carlo, and their branching fractions are used with the efficiencies to 
determine the size of each. Fortunately, the extra track and 7r° vetoes greatly reduce many 
backgrounds, such as tt+tt+tt^. Overall, we find a signal yield of about 2000 events. 

Although Figure HI shows all tag modes together, we actually fit each tag mode separately 
and calculate a branching fraction for each. The 6 branching fractions are then averaged to 
produce the final result. Table |V| shows the yields and efficiencies for each tag mode and 
the resulting branching fractions, without any systematic uncertainties or corrections. 

The systematic uncertainties are listed in Table IVII A small correction is applied for the 
particle identification of the vr"*" in — > i^'^vr"'", in addition to the uncertainty. The largest 
systematics arise from whether we allow the signal peak width to vary, from the extra track 
and extra vr" vetoes, from the shape of the signal peaks, and from the statistical uncertainty 
and input branching fraction of the Kgii^ background. The veto systematics include 

uncertainties on finding real tracks and vr^'s in background events, as well as on finding 
fake tracks and tt^'s in signal events. The fake track and 7r° systematics are determined by 
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Tag mode 


Tag Efficiency 
Factor 


Efficiency (%) 


D yield 


D+ ^ Kl'K+ 
yield 


Branching 
fraction (%) 


D- K+Tr'TT- 


0.9949 ± 0.0032 


82.31 ± 0.19 


80108 ± 342 


967 ± 37 


1.459± 0.056 


D- K+TT-TT-TT^ 


0.9579 ± 0.0055 


81.55 ± 0.27 


24391 ± 315 


345 ± 22 


1.662± 0.108 




0.9908 ± 0.0039 


82.37 ± 0.21 


11450 ± 144 


132 ± 14 


1.387± 0.147 


D- K^TT-TT^ 


0.9565 ± 0.0062 


81.94 ± 0.37 


25494 ± 404 


323 ± 23 


1.479± 0.108 


D- K^TT-TT-TT+ 


0.9552 ± 0.0050 


81.33 ± 0.24 


16739 ± 314 


184 ± 16 


1.291± 0.114 


D- K+K'TT- 


0.9870 it 0.0038 


81.09 ± 0.21 


6892 ± 154 


72 ± 11 


1.271± 0.195 


Sum or average 




81.80 ± 0.09 


165074 ± 723 


2023 ± 54 


1.456± 0.040 



TABLE V: Results for K^7r~^. The branching fraction for each tag mode is calculated 

from the corresponding yields, efficiency, and tag efficiency factor. This table does not include 
systematic uncertainties or corrections. 



Pion tracking 


± 0.35% 


Pion particle ID 


0.30 ± 0.25% 


Tag reconstruction: signal vs. non-signal 


± 0.2% 


vs D~ tags 


± 0.5% 


Kg veto systematics 


±1.1% 


Peak shapes 


±0.69% 


Fake D~ background shape 


±0.15% 


Fake D~ background yield 


±0.35% 


Background yields (except K^tt^) 


±0.49% 


— KgTT'^ efficiency Sz statistics 


±0.80% 


Fixed vs. floating peak width 


±1.63% 


Tail of signal peak 


±0.25% 


Total 


±2.43% 


— > KgTr~^ branching fraction 


±0.62% 



TABLE VI: Systematics for B{D+ i^£7r+). The "total systematic" does not include the D+ 
Kg7r~^ branching fraction systematic. 

looking for extra particles in fully-reconstructed DD events, in both data and Monte Carlo. 

The branching fraction, with systematics, is B{D+ Kl7r+) = (1.460 ± 0.040 ± 0.035 
± 0.009)%. The final uncertainty is the systematic uncertainty due to the input value of 

B. L>o ^ i^£7rO 

For the tagged D° — > K^n^ branching fraction measurement, the same 3 decay modes 
were selected with the same requirements as in the tagged — > Kgn'^ study. We require 
that there are no tracks and only one desired 7r°, and no rj on the other side. The invariant 
mass of the 7r° was required to be within 4 standard deviations of the known 7r° mass (same 
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FIG. 2: Missing mass distribution in data after removing events with extra tracks, vr'^'s, or r/'s. 

as used before). The invariant mass of rj was required to be within 3 standard deviations of 
the known t] mass. After rejecting the events with any rj or track or more than one vr*^, we 
compute the missing mass squared using the momentum of the D and 7r°, with both the D 
and 7r° masses constrained. In Fig. |21 we present the missing mass plots in data. 

Since D° K^tt^ only has one observable neutral particle, 7r°, the branching fraction 
bias for K^n^ is more apparent than that for K^n^. A correction factor was 

applied when computing the branching fraction. 

A number of background channels appear in the missing mass squared plot - K^tt^, 
r/vrO, tx^ti^, D° K*tt^, and "the rest" - with the tx^ti^ peak on the left 

side of our signal peak, and the K^^n^ and rjTv^ peaks right under our signal 

peak, as shown in Fig. |21 The total backgound is about 10% in the signal region. 




FIG. 3: Background for L>" K^vr". Upper left is K^tt", upper middle is 7]tt , upper right is 
lower left is K*Tr^, lower middle is "the rest", and lower right is all cases stacked together. 

In order to get the signal and estimate the background, we define three regions in M'^^^^-: 
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p-region (-0.1 ~ O.lGeV^), s-region (0.1 ~ O.SGeV^), and b-region (0.8 ~ 1.2GeV2). The 
backgrounds were split into three groups: tt'^tt"; Kgir^ and — > rjiT^; and "the 

rest". For K^tt^ and r^vr^, we have no experimental handles. We just trust the 

Monte Carlo, and use it for the subtraction. We use the yield in the p-region to estimate 
the background from tt^vt^, and the yield in the b-region to estimate "the rest." 

TABLE VII: Branching fraction for K^ir^. "s-s" means sideband-subtracted for tags, and 

background-subtracted for signal yield, "a" is the correction factor for tag bias. Error is statistical 
only. The Monte Carlo input branching fraction for — > K^tt^ is 1.06%. 





Ta 


gs 


Signal Yield 


MC 


Corrections 




MC 


raw 


s-s 


raw 


s-s 


eff(%) 


a 


BR(%) 


Ktt 


791054 


787214 


5415 


4907.5 


57.97 


1.00 


1.075±0.017 


Kmr^ 


1102536 


1066829 


7378 


6554 


55.36 


1.037 


1.070±0.015 


K3tt 


1171872 1136140 


7390.5 


6628 


52.37 


1.057 


1.054±0.014 


Data 


raw 


s-s 


raw 


s-s 


eff(%) 


a 


BR(%) 


Ktt 


48095 


47440 


367 


334.8 


57.97 


1.00 


1.217±0.073 




68000 


64280 


414.5 


363.1 


55.36 


1.037 


0.984±0.058 


KStt 


75113 


71040 


466.5 


418.0 


52.37 


1.057 


1.063±0.058 



After subtracting all the backgrounds, we get all the yields and computed branching 
fractions in Table IVTTl 

Contributions to systematic uncertainty are categorized in Table IVIIII The largest sys- 
tematic uncertainty comes from the extra 7r° veto. The difference of peak width and position 
between data and Monte Carlo produced the "Peak shape" systematic uncertainty. 



TABLE VIII: Systematic uncertainty (all in percentage) 



Systematic 


L>° ^ K-TT+ 


Z)0 ^ k-tt+tt'^ 


Z)° ^ K-TT+TT+TT- 


AE sideband 


0.14 


0.58 


0.57 


Background channel 


0.72 


0.98 


0.80 


Track simulation 


0.40 


0.0 


0.61 


Tag bias 


0.1 


0.1 


0.3 


Peak shape 


-0.90 


2.44 


0.75 


Extra vr'^ veto 


1.66 


1.66 


1.66 


7] veto 


-0.33 


0.33 


0.75 


vr*^ efficiency 








All (without TT^ efficiency) 


2.09 


3.18 


2.30 



We have measured the braching fraction of K^tt^ in the tagged data sample, and the 
branching fraction of K^ir'^ in both tagged and untagged samples. We use the two branching 
fractions of K'^ir^ to get the correction factor, 1 — 2rfCOs6 + r'j, and then apply the correction 
factor, 1 + 2rfCos6 + rj:, to the branching fraction of K^ix^ to get the true branching 

fraction of K^n^. 

rj ~ Rws is taken from recent Belle results 0, 0. By using the PDG value of "y", the 
branching fraction for — > i^^7r° was corrected to be 1.202 ± 0.016 ± 0.039%. Then we 
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calculate (1 — 2rcosS + r'^) for each tag mode separately. After correcting the branching frac- 
tions for the three tags separately and combining the results together, we get the branching 
fraction for D'^ ir£7r°: 0.940 ± 0.046 ± 0.032%. Note that the 7r° systematic uncertainty 
is not included. 



V. ASYMMETRIES BETWEEN D K'^tt AND D K^tt 

To compare D — > Kgir and D — > KItt, we compute the asymmetries 

B{D ^ Kin) - B{D ^ KItt) 
^ ' ~ B(D ^ K'^tt) + B{D ^ Kin) 

The error propagation in the asymmetry is complicated by correlations between the branch- 
ing fractions. For example, the D — > K^n measurements include an input D — > K^n 
branching fraction, so the two branching fractions are anti-cor related. Also, the mea- 
surements both include the same 7r° systematic, so this systematic cancels. 
The asymmetry is 

R(D+) = 0.030 ± 0.023 ± 0.025. 

The uncertainties are dominated by the — > Kgn~^ measurement, which used only a 56 
pb~^ subset of the 281 pb~^ data set. We expect an updated result soon, so the uncertainties 
on the asymmetry will improve significantly. 
The asymmetry is 

R{D^) = 0.122 ± 0.024 ± 0.030. 

This systematic uncertainty will also improve when the D^D^ cross section measurement is 
updated with the full 281 pb~^ data set. 



VI. INTERPRETATION 

The asymmetry measurements allow a measurement of the strong phase between A{D^ — > 
K~^n~) and A{D^ K~n^) under reasonable theoretical assumptions. 

The three Cabibbo-favored D Kn decays are described by an isospin 1/2 amplitude 
/1 1/2, an isospin 3/2 amplitude /I3/2, and their relative phase 5/. The amplitudes are 



^(D^ ^ K-n+) = J -Ay 2 + \ 0^3/26-*^^ (2) 




A{D' ^ irV) = -y + ^^^3/26-**^ (3) 

A{D+ ^ ^°7r+) = ^^3/26-^^^ (4) 

Without loss of generality, we may take A1/2 and ^43/2 to be real and positive. 

The four doubly-Cabibbo-suppressed D — > Kn decays are described by one isospin 3/2 
amplitude S3/2, and two isospin 1/2 amplitudes B1/2 and C1/2. It is reasonable to assume 



12 



that -B3/2 and are relatively real, and that S1/2, C'1/25 and A1/2 are relatively real. We 
can then write the amplitudes for the four doubly-suppressed decays 



AiD^ - 




A{D'- 




A{D+- 




A{D+- 





1/2 



c 



1/2J 



3/2^ 



2(^1/2 + C'1/2) 



'^i?3/2e- 



■1(5/ 



^(51/2 + 1/^53/26-^'^ 



(i?l/2 — Ci 



/2J 



3/2^ 



-1(5/ 



(5) 
(6) 
(7) 
(8) 



In our notation, A1/2, ^3/2, -B1/2, C1/2, and -B3/2 are all real, and the phase in any ampli- 
tude comes from 5/. We thus have six parameters to describe seven decays. A complete fit is 
underway, but results are not yet available. Instead, we illustrate what will be forthcoming 
by making some approximations. 

The exact expressions for the asymmetries (defined in equation ^ are 



R{D-^) 



2[Ai/2(5i/2 + C1/2) + 2A3/253/2 - V2(Ai/253/2 + ^3/2(51/2 + C1/2)) COS 5^ 



|Ai/2 - 72^3/26-^^^12 + |(5i/2 + C1/2) - V2S3/2e-*^^|2 
2^3/2(^3/2 + V^(5i/2 - C1/2) COS 5,) 



■(9) 



3^3/2 + 1 v 



1/2 



/2J 



1 R , ,= -i'5r|2 
3^3/26 ^1 



(10) 



From the three Cabibbo-favored decay rates, one readily shows that ^3/2 ~ (l/4)Ai/2 
and 6j ~ 90°. Furthermore, we assume that \A{D — >■ K^7r)\'^ is negligibly small compared to 
\A{D ^°7r)p. Thus, in the denominators of the expressions for R{D), the second terms 
(with only S's and C's) are ignored. We make these approximations in what follows. The 
asymmetries simplify to 



i?(D°) 



2(-Bi/2 + C'1/2) + -S3/2 
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8-B3/2 

3Ai/2 



A, 



/2 



(11) 

(12) 



The ratio of doubly suppressed to allowed decays r'j^^ = ^poZ^^-^+j is given by 



K-n 



1^1/2 



(13) 



and this simplifies to 



32 
33 



Bl/2 + C'1/2 



/2 



+ 



1 (B 



3/2 



^1/2 / 2 \A112 J 
Combining the (approximate) expressions for R{D^), R{D^), and r^^, we have 



.2 
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[3R{D 



_ 2i?(D°)i?(D+) + /?(D+)2] 
13 



(14) 



(15) 



^0.3 



+ 




-0.3 -0.2 -0.1 0.1 0.2 0.3 

Asym (D^) 

FIG. 4: Allowed values of R{D^) and R{D^). The size of the ellipse is proportional to rxn = 
y^B{D^ — > K~^-K^)/B{D'-^ K^TT^). Our measurements are shown, with uncertainties. 



Thus, and R{D^) must lie on an ellipse, whose size is set by r\^. 

Taking = 0.00363±0. 00038 from PDG 2004 branching fractions, we obtain the ellipse 
shown in Figured Our measured values for R{D^) and R{D~^) are shown, and they lie on 
the elhpse, as they should. 
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A. The D° K^TT^ Strong Phase 



We can also determine the 
of the two amphtudes 



K^Ti^ strong phase, defined as the phase of the ratio 



K+n-) 



re 



In particular, 



re 



Ostrong 



A{D^ K-7r+) 

_^(^l/2 + ^^l/2) + y|53/2e-* 



(16) 



(17) 



The simplified expression for (^strong is 



A 

'-'strong ~ 



tan 



1 
72 



/ 1 -5l/2+C'i/2 
4 ^1/2 



3/2 



M/2 



B1/2+C1/2 _^ 1 ^■ 



3/2 



Substituting :?1Z|±^ 
° ^1/2 



^i?(i50) - ^i?(D+) and 



^1/2 



8^1/2 



(18) 



3/2 



^1/2 



-"strong 



tan 



J_ ( R{D^)-3R{D^ 
71 V47?(i:'0) - R{D^ 



(19) 



Using the measured asymmetries as input, we find that the strong phase is consistent 
with zero: 



A 

^strong ~ 



(3 ± 6 ± 7)° 



The quoted uncertainties do not include uncertainties due to the approximations ^43/2 ~ 
(l/4)Ai/2 and Si 90°. These additional uncertainties would be at most a few degrees. 

We have performed a preliminary study of fitting for the amplitude parameters with the 
D — > Ktt rates as input. This fit produces results consistent with the above approximations 
on ^3/2/^1/2 and Sj. We note that we can relax one assumption by allowing one additional 
phase to be non-zero - for example, the phase of C1/2 relative to -B1/2. Studies of these more 
general fits are in progress. 



B. Constraint from D+ K+tt^ 

Consider the ratio of widths of the doubly-suppressed decays to charged kaon, I?+ — > 
K+tt'^ and ^ K+Tr': 

Naively, one would expect this ratio to be 1/2. Using the amplitudes given above. 



^ iyf(Sv2 + Cv2) + vl^3/2e-^^^P 

The naive result, p = 1/2, follows most simply from IC1/2/-B1/2I ^ 1, I-B3/2/-B1/2I ^ 1 - 
i.e., Ci/2 = -B3/2 = 0. We could also have I-B1/2/C1/2I <S 1 instead. A prehminary CLEO-c 
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result, B{D+ K+ti^) = (2.25 ± 0.36 ± 0.15 ± 0.07) x 10"^ Q, gives p = 0.64 ± 0.12, 
consistent with the naive expectation. 

Making the approximation 5/ ~ 90°, equation (PT|) simphfies to 



1 {Bi/2 — Ci/2y + 2Bl 
\2TTW 



^ 2 {By2 + Cy^Y + \Bl, ^^^^ 



We can rewrite this as 



^^ H^ -51/2+01/2^ ^^1/2+^1/2^ ^23) 

2 1 , 1 / -B3/2 



Call 



2 \ -Bl/2+C'l/2 



^^/^ - a (24) 



Bl/2 + (^1/2 
Bl/2 + C'l/2 



= /3 (25) 



Then equation (j23|) becomes 



1 (1 - 2l3f + 2a' 
" " 2 1 + la' 

Using the previously derived approximate expressions for R{D^) (jllj) and R{D^) (|12j). 
we find 

i?l/2 + Ci/2 3/?(D0) - i?(D+) ^ ^ 

Taking our measured values R{D+) = 0.030 ± 0.034 and R{D^) = 0.122 ± 0.038, we have 
a = 0.18 ± 0.24, which leads to = OmtoH 
Treating as a small number, we have 

2p ^ (1 - 2/3)2 + ^a' (28) 

Using the values of p and a, 

11-2/31 = 1.11^°:}^ (29) 

So P must be close to either or 1. If it is close to 0, /3 = — 0.06to;o5. 

If P is close to zero, then C1/2 <^ -Bi/2- If it is near one, then C1/2 ^ -Bi/2- One needs 
theoretical arguments to decide between these two cases. 
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